Possibility to disentangle anisotropic flow, flow fluctuation, and nonflow 
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We suggest the possibility to disentangle anisotropic flow, flow fluctuation, and nonflow using 
two-, four-, and six-particle azimuthal moments assuming Gaussian fluctuations. We show that 
such disentanglement is possible when the flow fluctuations are large, comparable to the average 
flow magnitude. When fluctuations are small, the disentanglement becomes diflicult. We verify our 
results with a toy-model Monte Carlo simulation. 
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I. INTRODUCTION 

Azimuthal distributions of hadrons produced in non- 
central heavy-ion collisions are anisotropic 1-4]. The 
distributions are often expressed by Fourier series, 



dN/dcj) cx 1 + 2vn cos(n</)) 



(1) 



n=l 



where cj) is particle azimuthal angle relative to the re- 
action plane. Coefficients w„ quantify the azimuthal 
anisotropics of different harmonics. Since reaction plane 
is not experimentally accessible, experiments use parti- 
cle azimuthal correlations to measure anisotropies @ . 
As a result, the measured anisotropies using two-particle 
correlation method, ^{2}, contain not only flow but also 
flow fluctuations and nonflow 

Flow fluctuations are mainly due to initial geometry 
fluctuations Q . Flow fluctuations could also result from 
fluctuating responses of the collision system to an identi- 
cal initial condition. Flow fluctuations cause a variation 
of anisotropy from event to event in a sample of sup- 
posedly the same events. We note, however, experimen- 
tally one cannot ensure the selection of the same events a 
priori, so flow fluctuations come also partially from real 
differences in the collisions, e.g. of different impact pa- 
rameters Q. 

Nonflow is caused by particle correlations that are 
not directly related to the reaction plane, such as jet- 
correlations and resonance decays 0] . In general, nonflow 
is a few-body correlation while flow is a many-body corre- 
lation. It was suggested that nonflow is largely cancelled 
in a particular combination of two- and four-particle cor- 
relations [13, [3- The resultant anisotropic flow from 
this four- particle method, w{4}, contains therefore sig- 
nificantly less nonflow than that from the two-particle 
method, v{2} 

Elliptic flow, V2, arising from the elliptical overlap re- 
gion of the colliding nuclei, has been extensively stud- 
ied Significant efforts have been invested to 
study V2 fluctuations and nonflow [l9l - [23| . Triangular 
flow, W3, can arise from the triangularity of the initial 
overlap geometry [23 - l30j , and may be dominated by fluc- 
tuations. So far flow fluctuations and nonflow have not 



been experimentally distinguished [23, [231 ■ 1^ well 
known that the two- and four-particle azimuthal mo- 
ments are not enough to determine the average flow, 
flow fluctuation, and nonflow [l^l- In this article we 
argue that it may be possible to determine the three 
quantities by the additional information of six-particle 
azimuthal moment. We support our argument by a toy- 
model Monte Carlo simulation. 



II. AZIMUTHAL MOMENTS 



The two-, four-, and six-particle n*'^ azimuthal mo- 



ments are jl2l |: 

(2) 
(4) 
(6) 



Jn(<t>i+(l>j-cl>k-4'i)^ 
Jn(4>i+4>j+cf>k~4'i-4>iJ.-4>v)\ 



(2) 
(3) 
(4) 



where the mean (...) is taken within a single event, and 
the indices must be all different i^j^/c^^y^/i^z/. 
The two-particle azimuthal moment is composed of flow 
and nonflow (S) [13, [11: 



(2) 



(5) 



The event-average two-particle azimuthal moment is the 
two-particle flow anisotropy measurement: 



v{2r ^ ((2)) = {v^) + {S) 



(6) 



To obtain the higher cumulants below, we assume flow 
and nonflow are uncorrelated and nonflow fluctuation is 
negligible. The four-particle azimuthal moment is (l3| 



(4) = + Av^S + 25^ . 
The four-particle cumulant flow is given by [l2j 
v{AY ^ 2((2))2 _ ((4)) = 2{vY - {v') . 
The six-particle azimuthal moment is [13, [13, HH 



(6) = 



9v*6 + ISv^S'^ + 66^ 



(7) 



(8) 
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The six-particle cumulant flow is given by [12|, |l3|, |3l| : 

v{6r ^ (((6))-9((2))((4)) + 12((2))3)/4 

= ((«6)-9(«2)(t;4)+12(^;2)3)/4_ (10) 

We note that the left sides of Eqs. ([S]) and pU)) are no- 
tations by definition because the right sides can be neg- 
ative. When the right sides are positive, then w{4} and 
v{6} are the four- and six-particle flow. 



III. GAUSSIAN FLUCTUATIONS 
ASSUMPTION 



Under the Gaussian fluctuation ansatz, we have 



(12) 
(13) 



-45(w)V^ + 15ct^ (14) 



We now assume that the fiow fiuctuations are Gaus- 
sian, 



dN 
dv 



1 



27rcr 



exp 



2f72 



(11) Substituting them into Eqs. 
I 



and ([TOl) . we obtain 



v{A}^^2{{2)f-m ^ {{vf-a^f~2a' 
i;{6rEE(((6))-9((2))((4)) + 12((2»3)/4 = {v)\{vf - ia^) . 
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(15) 
(16) 
(17) 



Eq. (ITC)) is as same as the result in Ref. [ifl, [III if the 
flow fluctuations are small so that ct* can be neglected. 
Under such circumstances, the flow fluctuation effect in 
ti{2}^ is -|-cr^ (positive) and in w{4}^ is — cr^ (negative). 



as seen from Eqs. ([T5|) and p6)) . respectively. 

In general flow fluctuation effects differ in different or- 
der azimuthal moments. Thus high order moments, com- 
ing with the unknown fluctuation effects, do not generally 
help resolving flow, flow fluctuation, and nonflow. When 
flow fluctuations are small, tT^/(w)^ ^ 1, even under the 
Gaussian fluctuation ansatz, higher order moments do 
not add new information. This is because Eqs. (jl6|) and 
([TTl) would give the same information: 

v{A}/{v) w H6}/(w) « 1 - aV2(w)' , ii<T^/{vf < 1 . 

(18) 

Only up to the order of / {v)^ do v{6} and v{4} differ: 



v{&}/v{A) « 1 - a^mvf , ii(T^/{vf < 1 



(19) 



Indeed data indicate that the six-particle W2{6| approxi- 
mately equals to the four-particle W2{4} p^.l32j. 

On the other hand, if flow fluctuations are large, then 
Eq. ([TTl) gives extra information than Eq. (fT6|). Figure [T] 
illustrates this point where the ratio of ti{6}/u{4} is plot- 
ted against cr^/(u)^: 



w{6}/w{4} 



(l-3aV(7;)2)i/6 



(20) 



As seen, when (t^/(u)^ ^ 0.2, a reasonable accuracy mea- 
surement of w{6}/v{4} would be able to determine the 



flow fluctuations. Note (T^/(w)^ is roughly the relative 
difference between w{2} and w{4}. 
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FIG. 1: Ratio of «{6}/u{4} as a function of /Iv)'^ under 
the Gaussian flow fluctuation ansatz. 



Odd harmonic anisotropics {vi and W3) are dominated 
by fluctuations. The fluctuations and the averages are 
likely on the same order. In fact, the average vi and 

with respect to the reaction plane are zero, and this 
may imply that their fluctuations roughly equal to their 
averages with respect to their own harmonic axes. It is, 
therefore, hopeful that the triangular (and directed) flow, 
flow fluctuation, and nonflow can be uniquely determined 
by the three equations of Eqs. (HH]), (HH), and dTT]). We 
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note, however, that in the case of small flow magnitude 
our neglect of the high order nonflow contributions may 
not be justified. 

For even harmonic anisotropics, such as the V2, the 
fluctuations are likely much smaller than the averages. 
In this case, the six-particle moment does not add much 
extra information beyond what is already in the four- 
particle moment. Thus, one will not be able to solve for 
flow, flow fluctuation, and nonflow from only two equa- 
tions given by the two- and four-particle moments. One 
can hope to overcome this by increasing the event statis- 
tics so that small deviation in i;{6}/w{4} from unity may 
be measurable. However, it is possible that even with 
infinite statistics, a measurement of u{6}/w{4} may not 
be enough to determine fiow, fiow fiuctuation, and non- 
fiow. This is because we have neglected those nonflow 
terms that are order (s) of magnitude smaller than the 
leading flow terms, but they may not be smaller than 
the high order terms in flow fluctuations. This depends 
on the details of nonflow, such as the multiplicities of 
multi-particle clusters, their anisotropics, and etc. 

It is worthwhile to note that the reaction plane is not 
pertinent to the two- and multi-particle correlation mea- 
surements of flow. In our Gaussian fluctuation ansatz, 
the only assumption is that the flow is a Gaussian dis- 
tribution about some mean {v)\ with respect to which 
plane the mean (u) is is not central to our method. The 
idea of Gaussian flow fluctuation is not new. The pos- 
sibility of extracting flow fluctuations under the Gaus- 
sian assumption via comparison of higher cumulants was 
studied earlier by Voloshin [33^. With relatively small 
fluctuation of V2 , the extraction of flow fluctuation is dif- 
ficult. In this work we point out that in the case of large 
fiuctuations (such as in U3), extracting fiow fiuctuation 
may be possible under the Gaussian fiuctuation ansatz. 

In real data, if fiow is not a Gaussian distribution 
with respect to any particular plane, then obviously our 
method will not work. For example, it was shown in 
Ref. [l0| that Gaussian distributions of the x and y ec- 
centricities in the reaction plane yield a Bessel-Gaussian 
distribution of the eccentricity magnitude in the partic- 
ipant plane. If elliptic flow is exactly proportional to 
the participant plane eccentricity then the harmonic mo- 
ments of the order of six and higher would give iden- 
tical information as in the fourth order harmonic mo- 
ment However, given an identical initial condition, 
the system response during collision evolution must also 
fluctuate. Such response fluctuations could smear the fi- 
nal flow distribution towards Gaussian. If, on the other 
hand, flow is indeed Gaussian with respect to a particu- 
lar plane, then the mean flow extracted from our method 
should be the (w) with respect to that plane. For elliptic 
flow, it is likely that this particular plane is the reaction 
plane, and the fluctuating participant plane gives rise to 
majority of the flow fluctuation. It is recently shown 
by Qiu and Heinz ^33] that the assumption of Bessel- 
Gaussian fluctuations works well for collisions of impact 
parameter 6 < 10 fm but breaks down in more peripheral 



collisions, while the hypothesis of Gaussian fluctuations 
works better than the Bessel-Gaussian assumption for pe- 
ripheral collisions, but breaks down for central collisions 
of 6 < 5 fm. 

Experimentally, one may turn the argument around. 
If one measured w{6} < w{4}, then it could mean that 
the Bessel-Gaussian fluctuation model is not completely 
correct. It does not necessarily mean that our Gaus- 
sian fluctuation assumption is correct, but it could pro- 
vide constraints to flow fluctuation models. If, on the 
other hand, one measured w{6} > w{4}, then neither 
the Bessel-Gaussian nor the Gaussian fluctuation model 
would be correct. One would have to devise new fluctu- 
ation model to accommodate experimental data. 

IV. TOY-MODEL SIMULATIONS 

We verify our results by a toy-model Monte Carlo sim- 
ulation. Below are the ingredients of the simulation. 

• For each event, we generate a V2 value according 
to the Gaussian distribution of Eq. (ITT|) with (W2) 
and (72 • We do the same for U3 with (113) and (T3. 
We allow negative V2 and of course negative U3 . We 
require |w„| < 0.5. 

• We set the reaction plane angle to zero. We gener- 
ate triangularity axis angle (V'3) randomly between 
and 2tt. We generate N — Nc\2 particles with 
azimuthal angle = Q-2tt according to the flow 
modulation of Eq. ([T}: 

dN/d(t) cx 1 + 2w2 cos 20 -I- 2u3 cos 3(0 - V'3) . (21) 

• ^ci2 of those N — Nc\2 particles each has a partner 
with identical azimuthal angle (so total number of 
particles is N , and iVci2 < N/2). These correlated 
pairs yield a nonflow correlation, which is simply 
given hy 6 = 2Nci2/N{N — 1). There are no clus- 
ters of more than two particles. The total particle 
multiplicity N and the number of directly corre- 
lated pairs iVci2 are fixed for all events; there are 
no multiplicity fluctuations or nonflow fluctuations 
in our simulation. 

We note that the analysis of V2 and are independent 
of each other, so one may simulate events with only V2 
or W3 (with corresponding fluctuations and nonflow) . We 
have simulated several cases, with both finite V2 and W3, 
only finite V2, or only finite V3. We obtained the same 
results. 

Table U lists the inputs and our results from several 
simulations. The values of the multiplicity iV in Cases 
(i) and (ii) correspond to real data of roughly 30-40% cen- 
trality of 200 GeV Au+Au coUisions at RHIC [H]. Those 
in Cases (iii) and (iv) correspond to real data of roughly 
the top 5% centrality of 200 GeV Au-f Au collisions at 
RHIC j3^ . The V2 values in Case (i) and (iii) correspond 
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approximately those measured experimentally in the re- 
spective centralities [l^. We used a relatively small 
fluctuation for Case (i) and large V2 fluctuation for Case 
(iii). Note cr|/(w2)^ ~ 0.2 implies an approximately 20% 
difference between v{2} and w{4}. We set the iVci2 value 
such that the input nonflow and flow fluctuation are on 
the same order, 6 ^ . 

We use the Q- vector method [Sj to calculate the two-, 
four-, and six-particle azimuthal moments from the sim- 
ulated events. The calculated moments are listed in Ta- 
blelH Also listed are the calculated v{2}, w{4}, and u{6} 
values for comparison. 

Using Eqs. p^ . ([TE)) . and pT|) . we can solve for flow, 
flow fluctuation, and nonflow. We can either do it numer- 
ically, or use the minimization method. In order to ob- 
tain the proper errors on the solutions, we use the latter. 
We treat the statistical errors on the moments as inde- 
pendent. We compute the normalized x^, the quadratic 
sum of the differences between the left and right sides of 
Eqs. (IT5t . P^ . and ([T7]). We minimize the to obtain 
the solutions of (w), a, and 5; the minimum x^ should in 
principle be zero as is the case in our solutions. The so- 
lutions are tabulated in Table ID As seen from the table, 
the resolved flow, flow fluctuation, and nonflow compare 
well to the inputs. We found, when the flow fluctuations 
are small, (T^/(w)^ <C 1, that the errors on (u) and a are 
strongly correlated with each other, and the error on 5 
is strongly anti-correlated. On the other hand, when the 
flow fluctuations are large, the errors on (u), cr, and 5 are 
mostly independent. 

We simulated 10^ events for each case. As seen, the 
errors on a and 5 are the largest for Case (i). This is 
because of the smaller jv^ in this case. As mentioned 
earlier, small jv^ makes the four- and six-particle mo- 
ments not much different. This results in a correlated 
errors on and that are not well controlled and may 
become large. This in turn gives a large error on b. 

V. DIFFERENTIAL FLOW 

In this section we discuss how to obtain differential 
flow, e.g. as a function of pt or particle species. For 



easy discussion we will use as an example, but our 
discussion is general for any type of differential flow. 

In order to obtain differential flow, one can form az- 
imuthal moments with one particle within the interested 
Pt region and the other particle(s) over the unrestricted 
Pt region. These azimuthal moments are identical to 
those in Eqs. ©, ([3]), and (|H), except that the flrst az- 
imuthal angle (^i is replaced by that of the particle within 
the interested pt region, 



(2') = (e^f-^^-^^)) , (22) 

(4') — (^gi-n(<t>'i+<Pj-<Pk-4>l)\j ^ (23) 

(6') = (^Q^^^'t''i+'t'i+'i'i'-'i'i-'^i^-'i'")'^ . (24) 
Similar to Eqs. ©, © and (O, we have 

(2') = v'v + 5' , (25) 

(4') « v'viv"^ +25) + 25'{v'' + 5), (26) 
(6') « v'v{v'^ + &v'^5 + Q5'^) + 

i5'{v^ + Av^S + 25''). (27) 



Suppose the anisotropic flow within the interested pt 
region is v' with Gaussian fluctuation cr'. We will assume 
that v' is completely correlated with v. That is, v' and v 
in the same event must satisfy: 

{v'-{v'))/a'^{v~{v))/a. (28) 

The probability of an event having v' and v is given by 
Eq. (dll)- Analogous to Eqs. (IT2t. (jTS]) . and (ITl). we have 

{v'v) = {v'){v)+a'a, (29) 
(«V) = {v'){v){{vY + + 

3(TV((«)2+a2) , (30) 
{v'v^) = {v'){v) {{v)^ + lQ{v)'^a^ + lba^) + 

b<j'a{{v)^ + &{v)^a'^ + i<7^) . (31) 

Similar to Eqs. (HI]), (dH) and (d?]), we have 



((2'» - 
2((2'))((2))-((4')) « 
(((6')) - 3((2'))((4)) - 6((2))((4')) + 12((2')) ((2))2) /4 « 



{v'v) + ,5' 

{v'){v)+a'a + 5' , (32) 
2{v'v){v'^) - {v'v^} 

{v'){v)i{vr-a')-a'a{{vr+a^), (33) 

- 3{v'v){v'^) - 6{v'^){v'v^) + 12{v'v){vY) /4 
{v)'^{{v'){vf - 2{v')a^ - {v)a'a) . (34) 



Here we have taken the nonflow correlation effect between one particle in the interested px region and another par- 
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TABLE I: Toy-model simulation to verify our method. Several cases are studied: (i) N — 400, A'ci2 = 72, (112) = 0.07, (72 = 
0.03,^2 = 9 X 10"''; (ii) iV = 400,iVci2 = 72, {vs) = 0,^3 = 0.03,^3 = 9 x lO""*; (iii) N = 1200, Afci2 = 288, (vz) = 0.03, (J2 = 
0.02, 52 = 4 X 10"''; (iv) iV = 1200, iVci2 = 288, (vs) = 0.03, aa = 0.03, 5^ ^ 4 x 10"**. We simulated 10^ events for each case. 



Case 


(i) 


(ii) 


(iii) 


(iv) 


((2)) 


(6.698 ± 0.003) X 10"^ 


(1.803 ±0.001) X 10"^ 


(1.7002 ± 0.0008) X 10"^ 


(2.199 ±0.001) X 10"^ 


((4)) 


(7.514 ± 0.007) X lO"'^ 


(7.29 ± 0.02) X 10"'' 


(5.849 ± 0.007) X 10"® 


(1.127 ±0.001) X 10"^ 


«6» 


(1.139 ±0.002) X lO"'' 


(4.79 ± 0.03) X 10"* 


(2.955 ± 0.008) X 10"* 


(8.95 ±0.02) X 10"* 


v{2} 


0.08184 ± 0.00002 


(0.04246 ± 0.00002) 


(0.04123 ± 0.00001) 


(0.04689 ± 0.00001) 


v{4r 


(0.0618 ±0.0001)'' 


(-7.9 ±0.2) X 10"'^ 


(-6.8 ±0.9) X 10"* 


(-1.59 ±0.02) X lO"** 


v{6r 


(0.0615 ± 0.0003)** 


(0.0 ±0.1) X 10"^ 


(-2.4 ±0.3) X 10"^° 


(-1.46 ±0.08) X 10"^ 


(v) 


0.070 ± 0.003 


0.007 ± 0.009 


0.0300 ± 0.0009 


0.0300 ± 0.0006 


a 


0.030 ± 0.004 


0.029 ± 0.002 


0.0200 ± 0.0005 


0.0300 ± 0.0001 


S 


(8.4 ±6.1) X 10"* 


(9.1 ±0.2) X 10"" 


(4.0 ±0.8) X 10"* 


(4.0 ± 0.4) X 10"" 



tide in the unrestricted px region to be S' . 

Knowing the (u), a, and S for particles in the unre- 
stricted pt region, we can readily solve for {v'), a' and 
(5' from Eqs. ([5^. ([55]). and ([M]). 

VI. CONCLUSIONS 

Anisotropic flow measurements have provided valuable 
information about the early stage of relativistic heavy- ion 
collisions. The contamination of nonflow, however, has 
hampered precise extraction of the hydrodynamic prop- 
erties of the collision system. Nonflow and flow fluctua- 
tions have not been successfully separated in experimen- 
tal measurements. We argue, if the flow fluctuations are 
Gaussian, that it may be possible to disentangle flow, 
flow fluctuation, and nonflow by simultaneous measure- 
ments of two-, four-, and six-particle azimuthal moments. 
We demonstrate that this possibility is real when the flow 
fluctuations are sizeable relative to the average flow mag- 
nitude. The disentanglement may be difficult when the 
relative flow fluctuations are small. We have verifled our 



conclusions by a toy-model simulation. 

To conclude, our study suggests that it may be pos- 
sible to measure elliptic and triangular flows, flow fluc- 
tuations, and nonflows uniquely in experiment. If, on 
the other hand, real data analysis does not yield reli- 
able results, our study offers additional insight-It may 
mean that either the Gaussian fluctuation assumption is 
not fully applicable and/or the flow fluctuations are too 
small. 
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